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Abstract: We study the spectrum of the confining strings in four- dimensional 
SU(A^) gauge theories. We compute, for the SU(4) and SU(6) gauge theories for- 
mulated on a lattice, the string tensions au related to sources with Zat charge k, 
using Monte Carlo simulations. Our results are consistent with the sine formula 
o"fc/o" = sinfc-^/ sin for the ratio between au and the standard string tension a. 
For the SU(4) and SU(6) cases the accuracy is approximately 1% and 2%, respec- 
tively. The sine formula is known to emerge in various realizations of supersymmetric 
SU(A^) gauge theories. On the other hand, our results show deviations from Casimir 
scaling. We also discuss an analogous behavior exhibited by two-dimensional SU(A^) 
X SU(A^) chiral models. 
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1. Introduction 

Quantum Chromodynamics is a nonabelian gauge theory based on the gauge group 
SU(3). The mechanisms underlying many of its fundamental properties, such as 
confinement, chiral symmetry, topological effects and the axial anomaly, are under 
active investigation; they are being studied by different approaches, including nu- 
merical simulations of the theory formulated on the lattice, several models of the 
vacuum, as well as some recent proposals derived from M-theory and AdS/CFT. 
Many features of QCD can be better understood by extending the study to SU(A^) 
gauge theories with A^ larger than three. In particular the large-A^ limit, which is 
obtained keeping g^N fixed {g is the gauge coupling) [[H, is of considerable inter- 
est from a phenomenological point of view, and is one of our best nonperturbative 
means of investigating QCD (see e.g. the reviews 0, ^, ^ [^ and references therein). 
Indeed, this limit is expected to preserve qualitatively most nonperturbative features 
of QCD. 

Four-dimensional non-Abelian gauge theories exhibit confinement, i.e. static 
sources in the fundamental representation develop a linear potential characterized by 
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a string tension a. As pointed out in many studies, it is important to investigate the 
behavior of the system in the presence of static sources in representations higher than 
the fundamental one. This may provide useful hints on the mechanism responsible 
for confinement, helping to identify the most appropriate models of the QCD vacuum 
and to select among the various confinement hypotheses. 

SU(A^) gauge theories confine by means of chromoelectric flux tubes carrying 
charge in the center Z^r of the gauge group. A chromoelectric source of charge 
k with respect to is confined by a fc-string with string tension cr^ (o"i = o" is 
the string tension related to the fundamental representation). If < ka, then 
a string with charge k is stable against decay to k strings of unit charge. Charge 
conjugation implies ak = CN-k- Therefore SU{3) has only one independent string 
tension determining the large distance behavior of the potential for k ^ 0. One must 
consider larger values of N to search for distinct fc-strings. The spectrum of the 
fc-string tensions is then determined by the ratios 

R{k,N) = —. (1.1) 
cr 

It has been noted that stable /c-strings are related to the totally antisymmetric 
representations of rank k, and that in various realizations of supersymmetric SU(A^) 
gauge theories R{k, N) satisfies the sine formula R{k, N) = S{k, N) where 

S(*.iV).!5«). (1.2) 

sm(7r/A) 

R{k, N) has been computed for the Af = 2 supersymmetric SU(A^) gauge theory 
softly broken to A/" = 1 |T^, obtaining Eq. (|1.2|). The same result has been found 



also in the context of M-theory, and extended to the case of large breaking of the A/" = 
2 supersymmetric theory [|l^. The same formula has been recently rederived [|ll 



using a different setup, i.e. gauge/string duality, suggesting that in the Af = 1 
supersymmetric gauge theories the sine formula may be quite robust. The interesting 
question is whether the sine formula holds also in nonsupersymmetric SU(A^) gauge 
theories. The M-theory approach to nonsupersymmetric QCD, although it is still 
at a rather speculative stage, suggests that this may be so |T^, 0. However, as 
discussed in Refs. |T^, |^, corrections from various sources cannot be excluded, so 
that this prediction cannot be considered robust. 

As pointed out in Ref . |]^ , it is interesting to compare the fc-string tension ratios 
in different theories. The idea is that such ratios may reveal a universal behavior 
within a large class of models characterized by SU(A^) symmetry, such as SU(A^) 
gauge theories and their supersymmetric extensions. Therefore, according to this 
universality hypothesis, the /c-string tension ratios in four-dimensional SU(A^) gauge 
theories should be given by the sine formula ( p..2| ). This notion of universality for 
the behavior of the A;-string tensions might complement the one conjectured for the 
type of effective string theory describing confining strings in gauge theories 0, |^ . 
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Another interesting and suggestive hypothesis is that the fc-string tension ratio 
satisfies the so-called Casimir scaling law [Q, i.e. R{k,N) = C{k,N) where 



C(.,iV).^ff^ (1.3) 

is the ratio between the values of the quadratic Casimir operators in the rank-A; anti- 
symmetric and in the fundamental representations. Casimir scaling is satisfied on the 
one hand by the small- distance behavior of the potential between two static charges 
in different representations, as shown by perturbation theory up to two loops [[17|] , 
and on the other hand by the strong-coupling limit of the lattice Hamiltonian formu- 
lation of SU(A^) gauge theories [15], |16]. Interest in Casimir scaling was recently 
revived [|18|, |1^, it has been triggered by numerical studies of SU(3) lattice 



gauge theory |2^, |23[, which indicate that Monte Carlo data for the potential be- 
tween charges in different representations are consistent with Casimir scaling up to 
a relatively large distance, r ^ Ifm. 

The Casimir scaling law holds exactly in two-dimensional QCD. In higher di- 
mensions no strong arguments exist in favor of a mechanism preserving Casimir 
scaling from small distance (essentially perturbative, characterized by a Coulombic 
potential) to large distance (characterized by a string tension for sources carrying 
Zn charge); nor across the roughening transition, from strong to weak coupling. We 
will show explicitly that Casimir scaling does not survive the next-to-leading order 
calculation of the ratios R{k, N) in the strong-coupling lattice Hamiltonian approach. 

It is worth mentioning another simple model for the spectrum of confining strings: 
if the interaction between fundamental flux tubes were so weak that no bound string 
states existed, then the spectrum would be given by 

F{k,N) = Mm[k,N - k]. (1.4) 

Note that all the above hypotheses considered have the same large- limit, i.e. 

S{k, oo) =C{k, oo) = F{k, oo) = k, (1.5) 

which is the expected result, since no bound states should exist for = oo. Note 
also that 

S{k,N) = k + 0{l/N^) . (1.6) 

In this respect the sine formula is peculiar because there are no a priori reasons for 
the large- A^ expansion of the /c-string tension ratio to be even in 1/A^. 

Of course, it is possible, and even likely given the current state of the theoretical 
knowledge, that none of the above hypotheses is correct. Nevertheless, we believe 
that a study able to discard some of them and determining the size of the correspond- 
ing corrections would be already important for the understanding of confinement in 
SU(A^) gauge theories. 
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The issue of the /c-strings can be investigated numerically using the lattice formu- 
lation of SU(A^) gauge theories. Recent numerical results for R{2,N), obtained for 
= 4, 5 |19|, 0, show that R{2, N) <2; thus, as < 2a, indicating that flux tubes 
attract each other, and definitely discarding the hypothesis ( |1.4|) of free strings. The 
available estimates of R{2, N) are substantially consistent with both the sine and 
Casimir formulas, thus they do not allow one to exclude any of the two hypotheses. 
This is also due to the fact that the two predictions for k = 2 are numerically close, so 
that high accuracy is necessary to distinguish them. In particular, the most precise 



result for the ratio R{2, 4), reported in Ref.|jT9|, lies between the predictions from the 



sine formula and the Casimir scaling, and is consistent with both within two error 
bars. 

The aim of this paper is to further investigate this issue. We present results from 
Monte Carlo simulations of the SU(A^) lattice gauge theories with A^=4, 6 using the 
Wilson formulation. For = 4 two independent fc-strings are expected, including 
the fundamental one. For = 6 there are three. We anticipate here our final results 
for the fc-string tension ratios: 

i?(2, 4) = 1.403 ±0.015, (1.7) 
i?(2,6) = 1.72±0.03, (1.8) 
i?(3,6) = 1.99±0.07. (1.9) 

Moreover, we found no evidence for stable string states associated with the symmetric 
rank-2 representation, in accordance with general arguments. 

Figure |I] summarizes our results, comparing our MC results with the above- 
mentioned hypotheses of spectrum. We claim that SU(4) and SU(6) results show 
substantial agreement with the sine formula, and therefore with the universality 
conjecture for the spectrum of the confining strings in asymptotically free theories 
with SU(A^) symmetry. The sine formula (0) predicts 5(2,4) = V2 = 1.414..., 
5'(2,6) = 1.732..., and 5(3,6) = 2. Moreover, the results show deviations from 
a strict Casimir scaling, whose predictions are C(2,4) = 4/3, C(2,6) = 8/5 and 
C(3,6) = 9/5. 

Considering our results all together, we can state that the sine formula is con- 
sistent within an accuracy of approximately 1%. This fact should be relevant for the 



recent debate on confinement models, such as those discussed in Refs. 18, HIS, 20 



21] , p5| , [26| , |27| , Pq , |29| , pOj |31| , p2| . Of course, our numerical results cannot prove that 
the sine formula holds exactly, but they place a very stringent bound on the size of 
the possible corrections. At the same time, our results appear rather conclusive on 
the existence of deviations from the Casimir scaling. Casimir scaling may still be 
considered as a reasonable approximation, since the largest deviation is about 10% 
for i?(3,6). 
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Figure 1: Comparison of tiie various hypotheses for the /c-string ratios with the Monte 
Carlo results. 



Finally, it is interesting to note that the sine formula ( |1.2| ) also emerges in the 
context of the two-dimensional SU(A^)x SU(A^) chiral models (see e.g. Ref. as 
a general reference), ci-dimensional chiral models and 2(i-dimensional lattice gauge 
theories present interesting analogies. In particular, the relation is exact for d = 1, 
and one can prove that Casimir scaling holds for the masses of the bound states. In 
analogy with four- dimensional SU(iV) gauge theories, in two-dimensional SU(A^)x 
SU(A^) chiral models the Casimir scaling law holds for the small- distance behavior 
of the correlation functions related to different representations. Moreover, it also 
holds for the strong-coupling limit of the corresponding lattice Hamiltonians, but it 
is not satisfied at next-to-leading order for a generic choice of the lattice Hamiltonian, 



such as the one of Ref. [33|. On the other hand, the exact S-matrix, derived using 



essentially the Bethe Ansatz [0, shows that bound states exist only for the 
rank-fc antisymmetric representations, and the ratio of their masses satisfies the sine 
formula. 

The paper is organized as follows. In Sec. |^ we present the results of the Monte 
Carlo simulations of the SU(4) and SU(6) lattice gauge theories. The rest of the 
paper presents analytical results that help in providing a more detailed picture of 
the characteristic features of the potential between static charges in higher-rank rep- 
resentations. In Sec. we describe the computation of the strong-coupling expansion 
in the lattice Hamiltonian approach of SU(A^) gauge theories, to the first nontrivial 
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next-to-leading order. We show explicitly that Casimir scaling is violated by the 
corrections to the leading order. In Sec. ^ we discuss the analogies between chiral 
models and lattice gauge theories. The appendices are dedicated to a number of 
issues related to the Monte Carlo simulations of SU(A^) gauge theory at large N, 
such as the matching of the lattice couplings in the large- iV limit, the bulk transition 
observed at finite bare coupling for sufficiently large, and the severe form of crit- 
ical slowing down which characterizes the Monte Carlo dynamics of the topological 
quantities such as the topological charge, and which appears to follow an exponential 
law rather than a power law. 

Short reports containing essentially our Monte Carlo results for the SU(6) gauge 



theory, and some preliminary results for SU(4), have already appeared in Refs. ^6 



2. /c-strings in four-dimensional SU(4) and SU(6) gauge the- 
ories 

In order to investigate the behavior of the /c-string tensions in gauge theories, we 
performed numerical Monte Carlo simulations of the four- dimensional lattice SU(4) 
and SU(6) gauge theories using their Wilson formulation 

W = -Nf3 Yl Tr [U^{x)U,{x + ^i)Ul{x + u)Ul{x) + h.c] . (2.1) 



In our simulations we employed the Cabibbo-Marinari algorithm [^] to upgrade 
SU(A^) matrices by updating their SU(2) subgroups (we selected 6 and 15 subgroups 
respectively for the SU(4) and SU(6) cases). This was done by alternating micro- 
canonical over-relaxation and heat bath steps, typically in a 4:1 ratio. In the following 
we consider a sweep as the upgrading of all links of the lattice independently of the 
algorithm; thus a over-relaxation and heat-bath cycle takes 5 sweeps. In Tables |I| 
and 1^ we present some information on our Monte Carlo runs for the SU(4) and SU(6) 
cases respectively. We provide the coupling value 

7-5^ (2.2) 

(this rescaled coupling is more natural for large N , due to the fact that the large- 
limit of the lattice theory is obtained keeping 7 fixed), the size of the lattices and the 
number of measurements, reported as the ratio between the number of sweeps A^^s^ 
and the interval between two measurements N.^. Moreover we report the values of the 
mean field ^ and cactus |^ improved couplings, 7mf and 7cactus respectively. 



and the value of the lattice spacing in units of the square root of the string tension. 
As discussed in App. |^, these quantities are useful to compare lattice results at 
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0.2699(23) 
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0.24523 


0.1906 
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3.17 
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0.24850 


0.1947 


123 


X 24 


4308k/20 


0.2368(6) 


2.84 





344 


0.25176 


0.1987 


123 


X 24 


2018k/20 


0.2122(8) 


2.55 










163 


X 32 


3615k/20 


0.2160(8) 


3.46 
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X 32 


4674k/20 


0.1981(5) 


3.17 



Table 1: Data sets available for SU(4). 
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lattice 
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a^fo 




0.342 


0.24234 


0.1843 


83 X 16 


213k/10 


0.3151(6) 


2.52 








123 X 24 


520k/20 


0.3239(8) 


3.89 


0.344 


0.24455 


0.1875 


123 X 24 


727k/20 


0.2973(5) 


3.57 


0.348 


0.24897 


0.1935 


103 X 20 


592k/20 


0.2534(6) 


2.53 








123 X 24 


712k/20 


0.2535(6) 


3.04 


0.350 


0.25117 


0.1963 


123 X 24 


442k/20 


0.2380(6) 


2.86 


0.354 


0.25556 


0.2017 


123 X 24 


270k/20 


0.2103(5) 


2.52 



Table 2: Data sets available for SU(6). 



different values of A^. In total, the whole study took about 10 years of CPU on a 
Pentium III IGhz cluster. 

The couplings were chosen to lie in the weak-coupling region. This is important 
because the Wilson lattice formulations of SU(A^) gauge theories undergo a first order 
phase transition for N sufficiently large, as argued using various approaches, such 



as Monte Carlo simulations H3, H3l HI, 45], mean field calculations [|q, H^, reduced 



models We found evidence for a first order phase transition in the SU(6) case, 
for 7c = 0.3389(4). This issue is discussed in App. |B|, where some additional Monte 
Carlo results are presented. Therefore, in our simulations we considered 7 values 
larger than 7c. Moreover, in order to avoid getting trapped in unphysical metastable 
states, we always used cold configurations as the starting point of our simulations. 
On the other hand, in the SU(4) case the MC data of the specific heat do not show 
any evidence for a bulk transition, contrary to some expectations coming from mean 
field calculations WA and earlier Monte Carlo simulations E^, IHl. The crossover 



between the strong- and weak-coupling region is characterized by a pronounced peak 
of the specific heat at 7 ~ 0.325 (corresponding to /3 ~ 10.4), similarly to the SU(3) 
case where the absence of a bulk transition is well established. 

We used asymmetric lattices ijj' x T) with a larger time size, along which the 
wall-wall correlations of Polyakov loops were measured. For some values of 7 we 
performed simulations for two lattice sizes in order to check for finite size effects. 
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The lattice sizes L were chosen so that L^/a ^2.5, and for most of them L^/a ~ 3. 
This requirement ensures that finite size effects on /c-string ratios are neghgible, as 
can be seen by comparison of the results for different sizes (see also Refs. ||49|| ) . 

In our simulations we measured also the topological charge Q, by a cooling 



technique, see e.g. Ref. ||5^. A severe form of critical slowing down is observed in 
this case, which worsens with increasing A^. Estimates of the autocorrelation time 
tq for Q, obtained from a blocking analysis of the data, turns out to be consistent 
with an exponential increase: tq oc exp(c^CT) where ^o- = o"^^^^, with c ~ 1.7 for 
SU(4) and c ~ 2.5 for SU(6). As a consequence, the run for the largest value of 7 
we considered for SU(6), i.e. 7 = 0.354, did not correctly sample Q, presumably 
because it was not sufficiently long («i 300k sweeps). This dramatic effect was not 
observed in the correlators used to determine the fc-string tensions, suggesting an 
approximate decoupling between the topological and nontopological modes. Indeed 
a blocking analysis did not show significant autocorrelations in measurements taken 
every 10-20 sweeps. In App. y this issue is discussed in more detail. The critical 
slowing down shown by the topological quantities represents a severe limitation for 
numerical studies of the topological properties at large values of A^ using standard 
Monte Carlo algorithms. The results related to the topological properties will be 
reported elsewhere. 

The /c-string tensions are extracted from the large-time behavior of correlators of 
strings in the antisymmetric representations, closed through the periodic boundary 
conditions (see e.g. Refs. pT|, p!9|): 

Cr{t) = 5^(x.[P(0;0)] Xr[P{x,,X2^t)]), (2.3) 

where 

P{xi,X2] t) = n^3f/3(xi, a;2, X3; t). (2.4) 

[/(x; t) are the usual link variables, and Xr is the character of the representation r. 
In particular for the fundamental representation: 

Xf[P] = TTP, (2.5) 

for the antisymmetric representation of rank k = 2 

Xfc=2[P] = TrP2_(TrP)2, (2.6) 

and for the antisymmetric representation of rank k = 3 

Xk=3[P] = 2TrP3 -3TrP2 TrP + (TrP)^ (2.7) 

These correlators decay exponentially as exp(— rrifct) where is the mass of the 
lightest state in the corresponding representation. For a fc-loop of size L, the fc-string 
tension is obtained using the relation [ pl[| 

nik = (TkL - —. (2.8) 
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The last term in Eq. (|2.8|) is conjectured to be a universal correction, and it is 
related to the universal critical behavior of the flux excitations described by a free 
bosonic string Numerical results for the three- and four-dimensional SU{N) 
gauge theories with various values of A^, see e.g. Refs. [|19|, |5^, and for the 



three-dimensional Z2 gauge theory support a universal description of the flux 
excitations in terms of a free bosonic string. So it is reasonable to assume that this 
picture, and therefore the relation ( |2.8| ), be valid every time that a stable string 
state propagates, regardless of the gauge group or the representation considered. 
The comparison of the results obtained for different lattice sizes will give further 
support to this assumption. Note that Eq. (|2.8|) is expected to hold for sufficiently 



large values of L. Ref. |^ argues that a lattice size L satisfying L^/a > 3 should be 
sufficient to observe a behavior according to Eq. (p.8|) for the loop mass. (Although 
Ref. considered SU(A^) gauge theories with N = 4,5, the case N = 6 should be 



equivalent in this respect.) 

In order to improve the efficiency of the measurements we used smearing and 
blocking procedures (see e.g. Refs. |5^) to construct new operators with a better 
overlap with the lightest string state. The smearing procedure replaces every spatial 
link on the lattice according to: 

Ukix)^plukix)+as J2 U,ix)Ukix + j)U^jix + k)\ , (2.9) 
[ ±0¥fc) J 

where V indicates the projection onto SU{N) and the sum only runs on spatial 
directions. The blocking procedure replaces the spatial links with super-links Uk{x) 
defined on a lattice with lattice spacing 2a (except for the time direction) according 
to 

Uk{x) =V lUk{x)Uk{x + k) + af Uj{x)Uk{x+3)Uk{x+3 + k)U]{x + 2k) 

(2.10) 

The blocking procedure can then be iterated n times to produce super-links of length 
2"a. The coefficients as and «/ can be adjusted to optimise the efficiency of the 
procedure. We constructed new super-links using = «/ = 0.5, three-smearing, 
and a few blocking steps, according to the value of L, i.e. one for L = 10, two for 
L = 8, 12 and three for L = 16. These super-links were used to compute improved 
Polyakov lines. 

We used a standard blocking analysis to check for possible autocorrelations in 
the wall-wall correlators used to determine the masses. The masses were obtained 
from fitting the time behaviour of the folded 2-pt correlators: 

Ffoided(t) = \ {Fk{t) + Fk{Nt - t)) = A{e-^-' + e-"^'=(^*-*)). (2.11) 
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• k=1 - fitted string tension - 
■ k=2 - fitted string tension - 

J I I I I 

1 2 3 4 5 

^min 

Figure 2: The A;-string tensions determined by the fits to the corresponding wall- wall 
correlations, as functions of the lower bound of the fitrange, tmin, for SU(4) and 7 = 0.341. 
We also show our final estimates, which are represented by the the continuous lines. 

The statistical error on the fitted parameters was computed using a bootstrap pro- 
cedure. 

The choice of the fit range [tmin, ^max] is a dehcate issue. In a lattice computation 
of masses the choice of the fit range is a source of systematic error that is very 
difficult to control. On the one hand, the data at early times are still contaminated 
by heavier states. On the other hand, the masses are so large in lattice units that 
the relative error of data increases rapidly with the distance. A large value of tmin 
reduces the systematic error due to contamination of heavier states, but leads to an 
increase of the statistical error. A satisfactory compromise is reached when the two 
errors are comparable, or, more cautiously, when the systematic error is estimated 
to be negligible with respect to the statistical one. 

The high statistics we collected for the SU(4) lattice gauge theory allowed us 
to achieve a good control of the systematic error coming from the contamination 
of heavier states. As an example, we discuss in some detail the analysis of the 
data obtained for 7 = 0.341, using approximately 2 x 10^ measurements. In Fig. ^ 
we show the results for o and 02 as obtained varying tmin, together with the final 
estimates that we will quote later. The value of tmax is not critical in this respect, 
the data reported in the figure have been obtained using tmax = 8, but the results 
are essentially independent of tmax- We observe clearly a dependence on tmin- Our 
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Figure 3: Data for the wall-wall correlations whose large-distance exponential behaviors 
determine the /c-strings, for SU(4) and 7 = 0.341. 

estimates are taken when a plateau is reached, i.e. for tmin = 4 in the case of a, and 
^min = 3 for (72- This should ensure that the systematic error due to heavier states 
is at most comparable with the statistical one. In Fig. ^ we show the data for the 
wall-wall correlations corresponding to the fundamental and k = 2 antisymmetric 
representations, and the curves ( |2.11| ) with our best estimates of the parameters. 
Thus we believe that the results and the errors we quote for SU(4) should account 
for this systematic error. 

For the SU(6) lattice gauge theory we could not afford such a clean analysis 
because of the relatively limited statistics. In this case fits were typically performed 
in the range [2 — 4] for 7 > 0.348, and [1 — 4] for the other values of the coupling. The 
systematic errors were checked by comparing the outcomes of fits over different time 
ranges. In practice, we could not use data for distances larger than t = 5, which did 
not allow us to check the dependence on the fit range as satisfactorily as in the SU(4) 
case. Thus the values of the fc-string tensions may be still subject to a systematic 
error due to the contamination of heavier states. However, we note that the ratios 
R(k, N) turn out to be more stable with respect to the choice of the fit range (this is 
already apparent from Fig. We indeed found that the variations of the fit range 
yield consistent results within the statistical error. So the estimates and the errors 
that we finally report for the ratios R{k, N) should be reliable also in the SU(6) case. 
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Results for the fc-string tensions obtained in our simulations are reported in 
Tables |^ and ^ for SU(4) and SU(6), respectively. The relative errors on the ratios 
R{2, 4), i?(2, 6) and R{3, 6) are essentially determined by the uncertainty on a2 and 
(J3, since the estimates of a are in general much more precise. A bootstrap analysis 
on the data performed directly on the ratios R{k,N), once the fit range has been 
chosen, usually provides statistical errors that are smaller than the ones we report. 
However we do not consider them sufficiently reliable, since the systematic error due 
to the contamination of heavier states is controlled only within the statistical error 
of the /c-string tensions. The ratios R{2, 4), R{2, 6) and i?(3, 6) are plotted in Figs. H, 
1^, and P respectively, versus a^a to evidentiate possible scaling corrections, which 
are expected to be O(a^) apart from logarithms. To facilitate the comparison, in the 
figures we show the predictions of the sine formula (|1.2|) and of the Casimir scaling 
( |1.3| ). Before going into details, it is worthwhile to emphasize some common trends 
in the results. Data confirm that the finite-size effects on the fc-string tension ratios 
are small for L^/a > 2.5. In all cases the approach to scaling is satisfactory, and the 
pattern of the scaling corrections turns out to be similar in SU(4) and SU(6). 

Let us discuss in more detail the results for the SU(4) gauge theory. Comparing 
the results obtained at 7 = 0.344 for L = 12 and L = 16, we see that when L^/o' > 2.5 
the size effects are small, and within the statistical errors of our data for the fc-string 
ratios. The size effects are instead observable within our statistical errors on a and 
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lattice 




a'^(T2 


0-2/0- 






0.335 


12^ X 24 0.0876(8) 


0.1203(17) 


1.374(20) 






0.337 


16^ X 32 0.0728(12) 


0.103(3) 


1.413(38) 






0.338 


12^ X 24 0.0698(4) 


0.0973(12) 


1.395(18) 






0.341 


12^ X 24 0.0561(3) 


0.0786(10) 


1.402(17) 






0.344 


12^ X 24 0.0450(3) 


0.0634(7) 


1.407(15) 








16^ X 32 0.0466(3) 


0.0656(9) 


1.408(20) 






0.347 


16^ X 32 0.0392(2) 


0.0549(4) 


1.400(10) 








Table 3: fc-string tensions for SU(4). 
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lattice 






aVs 


0-2/0- 


0-3/0- 


0.342 


8^ X 20 


0.0993(4) 


0.164(1) 


0.190(3) 


1.65(2) 


1.91(3) 




12^ X 24 


0.1049(5) 


0.174(3) 


0.201(9) 


1.66(3) 


1.91(9) 


0.344 


12^ X 24 


0.0884(3) 


0.153(2) 


0.173(4) 


1.73(2) 


1.95(5) 


0.348 


10^ X 20 


0.0642(3) 


0.111(2) 


0.134(7) 


1.73(3) 


2.08(10) 




12^ X 24 


0.0642(3) 


0.110(3) 


0.132(7) 


1.71(4) 


2.06(11) 


0.350 


12^ X 24 


0.0567(3) 


0.097(2) 


0.110(6) 


1.72(3) 


1.95(9) 


0.354 


12^ X 24 


0.0442(2) 


0.0766(11 


) 0.090(3) 


1.73(3) 


2.04(6) 



Table 4: fc-string tensions for SU(6). 
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Figure 4: The scaling ratio i?(2,4) as a function of a^a for S'f7(4). 



(72. The data show a good scahng behavior. Only the results for the smallest value of 
7, i.e. 7 = 0.335 apparently show scaling corrections. We extract our final estimate 
for the ratio (T2/cr from the data at the largest values of 7, i.e. 7 = 0.341, 0.344, 0.347, 
and for their corresponding largest lattices. Combining them to obtain the central 
value and taking the typical error of each single point as estimate of the error, we 
propose as final estimate i?(2,4) = 1.403 ± 0.015. Of course, this estimate assumes 
that the scaling corrections are already small and negligible for a^cj ~ 0.05. The data 
for smaller 7, i.e. 7 = 0.335, 0.338, are used to check this hypothesis. They indicate 
that the scaling corrections are of the same size as the error we reported above. We 
will return to this point later. Our result is consistent with the prediction of the sine 
formula, S'(2,4) = = 1.414..., within an uncertainty of approximately 1%. On 
the other hand, this result does not support Casimir scaling, whose prediction in this 
case is C(2,4) = 4/3. 

i.e. i?(2,4) 



The result obtained in Ref. IS 



1.357(29), is marginally consistent 
with ours. The comparison improves considering the results i?(2,4) = 1.377(35), 
obtained in Ref. |]19| by discarding the data for the smallest value of /?, i.e. P = 10.55 
corresponding to 7 = 0.3297.... 

Let us now consider the results for the SU(6) gauge theory. Comparing the 
results for the string tension ratios for different lattice sizes at constant 7 shows 
little finite size effects. As in the SU(4) case, finite size effects are observed on the 
A;-string tensions, but they cancel out in the ratios. The ratio R{2, 6) displays good 
scaling for 7 > 0.342. Given the good scaling behaviour, again we do not attempt to 
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Figure 5: The scaling ratio R{2,6) as a function of a^a for SU{6). 

fit the dependence of our result on the lattice spacing a. Our final value for R{2, 6) 
is obtained averaging the results at 7 = 0.348 and 7 = 0.350, while the error is 
given by the typical error of each single point. Analogously to the SU(4) case, the 
data at smaller 7, and in particular the one at 7 = 0.344, are used to check that 
the scaling corrections are small, and at most comparable with the error we report. 
Given the poor sampling of the topological charge in the run at 7 = 0.354, we do 
not include this value of 7 in the final estimate of the string tension ratio. Similar 
comments apply to the R{3, 6) ratio. We finally note that the results for R{k, 6) at 
7 = 0.354 are in agreement with those obtained at smaller values of 7, for which 
Q was sampled correctly, lending further support to the previously mentioned large 
decoupling of the modes determining the string tensions and the topological ones. 
As anticipated in the introduction, our final estimates are R{2, 6) = 1.72 ± 0.03 and 
R{3, 6) = 1.99 ± 0.07. They are both consistent with the predictions of the sine 
formula ( [L.2[) , which are 5(2,6) = 1.732... and 5(3,6) = 2, respectively, within an 
uncertainty of approximately 2% for k = 2 and 4% for k = 3. On the other hand, 
they do not support Casimir scaling; the predictions in this case are C(2,6) = 8/5 
and C(3,6) = 9/5. 

We have also explored correlators in the symmetric rank-2 representation, whose 
character is given by 

Xfc=2,symm[P] = Tl + (Tl Pf , (2.12) 

finding no evidence for new stable string states for both SU(4) and SU(6) cases. The 
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Figure 6: The scaling ratio R{3, 6) as a function of a^a for SU{6). 



masses extracted in the symmetric channel were consistent with msymm > 2mf, as 
expected because the symmetric string should indeed decay into two fundamental 
strings. In this case, Eq. (|2.8| ) should not apply and we did not try to extract a 
string tension. 

In conclusion, our results turn out to be consistent with the sine formula, and 
show that there are sizeable corrections to the Casimir scaling prediction. In order 
to further check the robustness of our statement concerning the deviation from the 
Casimir formula, we have analyzed our data using also fits which allow explicitly 



for scaling corrections. As suggested in Ref. [0, one may fit the data to the linear 
behavior A + Ba. For this purpose one must consider all data, including the ones 
for the smallest values of 7, since they are the only ones showing apparent scaling 
violations. In the SU(4) case we find R{2, 4) = 1.424(23) with ^ 1-0. In the SU(6) 
case the results are i?(2, 6) = 1.76(4) (x^ ~ 2.6), /2(3, 6) = 2.13(11) (x^ ~ 1.4) and 
i?(2,6) = 1.77(6) (x^ ^ 2.7), i?(3, 6) = 2.16(19) (x^ ^ 1-3) respectively with and 
without the point at the largest value of 7, i.e. 7 = 0.354. Therefore, if one wants 
to be more cautious in treating the systematic error due to scaling corrections, one 
may take into account the difference between the linear fits and the nonextrapolated 
data, yielding 



i?(2,4) = 1.403 ±0.015 
i?(2,6) = 1.72±0.03 1°:°^, 



(2.13) 
(2.14) 
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R{3, 6) = 1.99 ± 0.07 toll, (2-15) 

which cover all the results obtained above. These conservative estimates are still 
not consistent with the Casimir formula. Thus, our conclusions are fully justified 
even by this overly cautious analysis. One should note that corrections to Casimir 
scaling are to be expected as discussed in the previous section (see also Sec. H). The 
implications of these results for models of the Yang-Mills vacuum deserve further 
investigation. 

One may write down a general expression for R{k, N) taking into account the 
following constraints: -R(l, A^) = 1 by definition, R{k, N) = R{N — k, N) by charge 
conjugation, and R{k, oo) = k which is the expected large- iV limit. A general ex- 
pression satisfying these conditions can be written as 



where 

t;(l-,JV) = sm — 



l + ^c,(iV) (si„^^)' 

i = l ^ 



(2.17) 



and Cj are coefficients which may depend on A^, but in such a way as not to spoil the 
large- limit, i.e. Ci{N)/N^ = 0{1/N). Our results show that, if the sine formula is 
not exactly satisfied, the corrections must be small, thus q <^ 1. In order to better 
quantify this statement, one may keep only the first term in the sum of Eq. ( p. 17] ) 
and find a bound on ci (assuming it constant). Our Monte Carlo results provide the 
bound |ci| < 0.05. We think that the accuracy of the sine formula in predicting the 
(Tk/c ratio should trigger further fundamental theoretical investigations. 



3. Strong-coupling expansion in the lattice Hamiltonian ap- 
proach of SU(A^) gauge theories 

The strong-coupling expansion in the lattice Hamiltonian formulation of the theory 
is an analytical approach that can be used to investigate Casimir scaling and its 
corrections. 

The lattice Hamiltonian of a /^-dimensional SU(A^) gauge theory is defined in 
terms of link operators in the fundamental representation ij^{x). Following Ref. 
we consider the Hamiltonian 



™ 2a 



where £"^(0 = Yl!l=i ^ E°-{l)E"-(l) is the quadratic Casimir operator of SU(A^) associ- 
ated with the link I of a. {D — l)-dimensional cubic lattice, and U{p) is the product of 
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link operators on the boundary of a plaquette p. The SU(A^) gauge theory is reahzed 
in the continuum hmit (7 — 0. The lattice Hamiltonian formulation is not unique, 
and therefore the corresponding strong-coupling expansion in powers oil/g should be 
considered as a suggestive investigation method, but it is difficult to extract reliable 
quantitative information from it. 

In the perturbative strong-coupling approach in which g 00, one works in the 
space of states Ylx 1^) diagonal in U: 

Uf\U) = Uf\U), (3.2) 
[E\Ur] = -\KU. (3.3) 

where r indicates the generic representation and A" the corresponding generators. In 
the strong-coupling limit g ^ 00, the vacuum |0) is the lowest eigenstate of 

Ho^Aj^E'il), (3.4) 

I 

where A = g'^/2a, thus it satisfies £'^|0) = 0. In order to compute the fc-string 
tensions, we must consider the force law between widely separated quarks in the 
strongly coupled limit. The potential energy is defined as the lowest energy compat- 
ible with the presence of a quark q at the origin and an antiquark q at site s. The 
minimum-energy gauge-invariant state is obtained by exciting the shortest path of 
links connecting the qq pair. Let us take the qq pair along a side of the cubic lattice, 
thus the least number of excited links is I = s. In the strong-coupling limit g 00 
it can be written as the product of I link operators joining the origin to the site s: 

|r,0 = d;'/^Ur{l)Ur{'^)...Ur{s)\0), (3.5) 

where r is the representation of SU(A^) considered and its dimension. The inner 
product (I) is defined using the invariant integration over the group: 

j dgUfUf,'' = d;'6rr'Saa'Sbb'- (3.6) 

The energy corresponding to the state |r, I) is given by the matrix element 
{r,l\Ho\r,l). As already mentioned, the leading order of the energy is proportional 
to Cr = |A"AJ?, the value of the quadratic Casimir operator in the representation r, 
i.e. 

{rJ\Ho\rJ) = AlCr. (3.7) 

Higher order corrections can be computed by systematic application of perturbation 
theory, where the perturbation interaction Hj is given by the second term of the 
Hamiltonian ( |3.1[ ). The perturbation Hj corrects the energy of these states and of 
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the vacuum to 0{1/ g^), i.e. to second order in perturbation theory. These corrections 
are obtained by computing 

1 



{r,l\Hi 



-Hj\rJ) 



(3.8) 



E — Hq 

and subtracting the corresponding second order contribution to the vacuum energy 

1 



{0\Hj—-Hj\0). 



(3.9) 



For the fc-string tension associated with the rank-/c antisymmetric representation 
(1*^; 0) (where l'^ indicates one column with k squares in the corresponding Young 
tableau) we obtain 



cTfc = AC, 



(lfe;0) 



l+(^e4(A:,iV) + 



1 



(3.10) 



where C(ife.o) is the value of the quadratic Casimir operator for the representation 
(1^0), and 



8N^ 

Cnk 



(lfe;0) 



dr 



2C(1;0) C'(ife.o) — 3C(i;0) — Cr^ Nd^ik.Q) 



(3.11) 



The sum extends over the four representations obtained by composing the rank-fc 
antisymmetric representation (1*^; 0) with the fundamental (1; 0) and antifundamen- 
tal (0; 1) representations, i.e. (l''+^;0), (2, l''"^; 0) and (l''"^;0), (l''; 1) respectively. 
For notation see e.g. Ref. [^. The corresponding dimensions and Casimir values 

Cf- are 



(i(ife.o) — 

C(2,l'''-i;0) 
dnk.^\ = d, 



N 
k 

k{N + 1){N - k) 



k 



2N 

N + 1' 
k + 1 

{k + 1)[N^ - N{k - 2) - k - 1] 



2N 



i(ife.i) - 



(2,lJV-fc-l;0), 



(2,1^ 



;0)- 



(3.12) 

(3.13) 
(3.14) 

(3.15) 
(3.16) 
(3.17) 



Notice that in the above expressions (l''; 0) = (0;0), i.e. the singlet representation. 
For k = 1, i.e. the fundamental string tension, the expression ( |3.10| ) reproduces the 
result reported in Ref. |5^, i.e. 



a = A- 



2N 



1 - 



D 



16N^{3N^ - 5) 



+ 



(3.18) 
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One can easily see that e^lk, N) satisfies the relation 64(^5 A^) = e^lN — k, N). It has 
the expected limit for N ^ 00, indeed the term coming from the subtraction of the 
vacuum contribution, which is the first one in Eq. ( |3.11| ), cancels the leading power 



in A^, and the expression goes to a constant for iV — > 00. 
For the fc-string tension ratio we obtain 

a,_kiN-k)f^ , [e,{k,N)-e,{l,N)]+0[^—-]\. (3.19) 



Explicitly for k = 2, 

64(2, iV) - 64(1, iV) = 

AN^{N - 3)(12A^6 + UN^ - 4QN^ - 47N^ + QIN^ + 34A^ - 24) 
{N - 2){2N - 3){2N + 3){N'^ - \)\2N'^ - 1)(2A^2 + _ 4) 

Note that the correction to the leading order Casimir scaling is 0(1/A^). 

This result provides an explicit example of corrections to the Casimir scaling 
prediction: while the leading order strong coupling calculation yields Casiir scaling, 
0(1/A^) corrections arise at the next-to- leading order. This fact holds for both D = 4 
and D = 3. 



4. Chiral models and lattice gauge theories 

In this section we discuss the analogies existing between 2(i-dimensional SU(A^) gauge 
theories and d-dimensional SU(iV)xSU(A^) chiral models. In particular, we compare 
the spectrum of the fc-strings in SU(A^) gauge theories with the known spectrum of 
the bound states in chiral models. 



4.1 Analogies 

Unitary matrix models defined on a lattice can be divided into two major groups, ac- 
cording to the geometric and algebraic properties of the dynamical variables: when 
the fields are defined in association with lattice sites, and the symmetry group is 
global, i.e., a single SU(A^)i:, x SU(A^)ij transformation is applied to all fields, we are 
considering a spin model (principal chiral model); in turn, when the dynamical vari- 
ables are defined on the links of the lattice and the symmetry is local, we are dealing 
with a gauge model (lattice gauge theory). An analogy between c?- dimensional chiral 
models and 2(i-dimensional gauge theories can be found according to the following 
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correspondence table |^: 

spin gauge 

site, link link, plaquette 

loop surface 

length area 

mass M string tension a 
two-point correlation Wilson loop 

To this correspondence table, one may also add: the spectrum of the bound states on 
the side of the chiral models and the spectrum of confining strings for SU(A^) gauge 
theories, i.e. 

spin gauge 

bound state mass fc-string tension ak 

While the above correspondence in arbitrary dimensions is by no means rigorous, 
there is some evidence supporting the analogy. 

Let us consider the following lattice formulation of (i-dimensional SU(A^) x SU(A^) 
principal chiral models 

^chiral = -2N(3 ReTr [U{x)U\x + /i)] (4.1) 

(where [3 = 1/NT, e SU(A^) and /i = l,...,c?); for SU(A^) gauge theories in 2d 
dimensions we consider the Wilson lattice formulation ( |2.1| ). 

In the case d = 1 one can prove an identity between the partition function (and 
appropriate correlation functions) of the two-dimensional lattice gauge theory and 
the corresponding quantities of the one-dimensional principal chiral model (see e.g. 
Ref. [§]). Both theories are exactly solvable and the correspondence can be explicitly 
shown. In particular, using the results of Ref. |Q, one may easily show that the 
Casimir formula holds for the spectrum. Thus 



Mk _ k{N - k) 
~W ~ N -I 

holds for the one- dimensional SU(A^) xSU(A^) chiral models, and 

ak _ k{N - k) 
V ~ N-1 



(4.2) 



(4.3) 



for two-dimensional SU(A^) gauge theories. 

Although for higher dimensions the correspondence does not strictly hold, we 
still have some analogies: 

(i) The two-dimensional chiral model and the four- dimensional non-abelian gauge 
theory share the property of asymptotic freedom and dynamical generation of a mass 
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scale. In both models these properties are absent in the Abelian case (XY model 
and U(l) gauge theory respectively). 

(ii) The classical equations of motion describing the dynamics of the spin vari- 
ables in the two-dimensional chiral model and the Wilson loops in the four- dimensio- 
nal nonabelian gauge theory are similar (see e.g. Ref. 0). It turns out that the gauge 
fields are chiral fields on the loop space. This analogy persists at the quantum level 
as well, since the classical equations of motion dictate relations among correlation 
functions in the corresponding quantum theory. 

(iii) Exploiting asymptotic freedom, one can use perturbation theory to de- 
termine the short- distance behavior of the potential of two heavy quarks in four- 
dimensional SU(A^) gauge theories and of the two-point correlation function in two- 
dimensional SU(A^)xSU(A^) chiral models. Both show Casimir scaling. In the four- 
dimensional SU(iV) gauge theory the force Fr{x) between two heavy quarks in the 
representation r is proportional to Cr up to two loops jl^. Resumming the lead- 
ing logarithms using standard renormalization-group arguments, the short-distance 
behavior of the force Fr(x) is: 



. X ^ ar{l/x) 

(X) = —Cr 



X^ 



ar{ylx) 



boL{x) 



L{x) \ L2 



where L{x) = — InxA^; and A^, is a mass scale (the so-called A-parameter associated 
with the above definition of running couplings ar{x); bo and bi are the first two uni- 
versal coefficients of the perturbative expansion of the /3-function (3{a) = —xdar/dx. 
bo, bi and A^,. are independent of Cr- A dependence on Cr of the O ((lnL)^/L^) term 
in ar{l/x) is not excluded. Its computation requires a three-loop calculation of the 
potential. As we shall see in Sec. ^4.2| , Casimir scaling emerges also in the short- 
distance behavior of two-point functions in chiral models. 

(iv) Another analogy concerns the large- limit, which in both cases is given by 
a sum of planar graphs. In this limit particles in chiral theories become free, indeed 



the known S-matrix ||3J, ^ becomes trivial. On the other hand, it is conjectured 
that the large- A^ limit of gauge theories is a free string theory, see e.g. P], although 
its nature is still not clear. 

(v) Their strong-coupling expansions look similar, provided that point-like exci- 
tations of chiral fields are substituted by closed fiux lines for the gauge fields. 

(vi) Approximate real-space renormalization recursion relations obtained by Mig- 
dal are identical for rf-dimensional chiral models and 2(i-dimensional gauge mod- 
els. 

4.2 The two-dimensional SU(A^) xSU(A^) chiral models 

The two-dimensional SU(A^)x SU(A^) principal chiral models are asymptotically free 
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matrix-valued field theories defined by the action 

S=^j d\TT d^U{x) d^U\x). (4.4) 

where U{x) G SU(A^). 

Using the existence of an infinite number of conservation laws and Bethe-Ansatz 
methods, the on-shell solution of the SU(A^) x SU(A^) chiral models has been proposed 
in terms of a factorized S-matrix |35|. The analysis of the corresponding bound 
states leads to the mass spectrum 

A4 = Af5l!^<^. l<k<N-l, (4.5) 

sm(7r/A') 

where Mk is the mass of the fc-particle bound state transforming as a totally antisym- 
metric tensor of rank k. M = Mi is the mass of the fundamental state determining 
the Euclidean long-distance exponential behavior of the two-point Green's function 

G(x) = l(Tr[/(OMx)t). (4.6) 

No bound states exist for other representations. Correlation functions associated 
with the generic representation r can be defined by 

Gr{x) = j^{Xr[um{xy]), (4.7) 

where dr and Xr are respectively the dimension and the character of the representa- 
tion r. The structure of the S-matrix implies that stable bound states propagate only 
in the totally antisymmetric representations with masses given by Eq. (|4.5| ). Note 
that, according to the correspondence discussed in Sec. ^A\ , the correlation functions 
Gr{x) play the role of the r-representation Wilson loops for SU(A^) gauge theories. 
The mass spectrum ( |4.5| ) has been verified numerically at = 6 by Monte Carlo 
simulations |^, |59| . 

As in SU(A^) gauge theories, the large- A^ limit of these models is formally rep- 
resented by a sum over planar graphs. The ^-matrix has a convergent expansion in 
powers of 1/A^, and becomes trivial, i.e., the ^-matrix of free particles, in the large- A^ 
limit. Note also that, as in the case of SU(A^) gauge theories, no a priori reasons 
exist for the large- A^ expansion of the fc-state mass ratios ([4.5|) to be even in 1/A^. 

Asymptotic freedom allows us to determine the small-distance behavior of the 
correlation functions in perturbation theory. In two-dimensional SU(A^) xSU(A^) 
chiral models, the logarithm of the two-point function Gr{x) is the analog of the po- 
tential of two separated quarks in the representation r. The small-distance behavior 
of In Gr{x) satisfies Casimir scaling, similarly to what happens for the potentials in 
four-dimensional SU(A^) gauge theories. Extending the results of Ref. pertaining 
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to the fundamental two-point function, in the MS renormahzation scheme associated 
with the X space we find 

N 



( iV 
hiGr{x) = —Crt< In(yux) + -—t\n{fix) 
I 8 

2C(3) 



(4.8) 



+t 



^(3 

256 V 



+ 



3N\ ^ , A^^ 

\n{fix) + — \n{fixy 



+ 



iV2 

48 



\n{fixf +0{t 



where Cr is the Casimir value of the representation r. As in four-dimensional SU(A^) 
gauge theories, one may define a running coupling tx from the relation 

d\nGr(x) „ 



. . (4.9) 

OX X 

Standard renormalizat ion-group arguments allow one to resum the leading loga- 
rithms, yielding 



1 



hL{x) 



hi In L{x) 
L(x) 



+ 



bi 



(4.10) 
A^V(1287r2) are 



As before, L{x) = — InxA^^ and A^; is a mass scale; bo = N/ (Sir) 
the first two coefficients of the /3-function Pit^) = —xdi^/dx = —botl — bit^ + 0(t^). 
We also mention that at the next order of perturbation theory, i.e. O(t^) in Eq. ( ^4.8| ), 
there are three-loop diagrams whose group factor would violate Casimir scaling. 
Thus, barring particular cancellations, Casimir scaling should not be satisfied by the 
O ((In L)^/L^) term of Eq. ( [4.10[ ). The same observation applies to four-dimensional 
SU(A^) gauge theories. 

Finally, again similarly to QCD, one may easily check that, using the lattice 
Hamiltonian approach of Ref . ||33| , Casimir scaling is recovered in the strong-coupling 
limit; it is violated by higher order corrections. 

The analogy with chiral models highlights some general trends: Casimir scaling 
is exact in d = 1, as it is in the case of SU(A^) Yang Mills theories in d = 2; for 
d = 2, both perturbation theory and strong coupling yield Casimir scaling to lowest 
order. However next-to-leading order calculations explicitly show corrections to such 
behavior. This fact is consistent with the picture that emerges for d = 4 SU(A^) 
gauge theories from our Monte Carlo simulations. 
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A. Rescaled and effective lattice couplings 



In order to compare Monte Carlo results for various values of A^, it is useful to 
introduce the rescaled coupling 



7 



glN 2iV2 



(A.l) 



which is kept fixed in the large- limit of the lattice theory, see e.g. Refs. 0, §]. As 
already noted in Refs. [0, the correspondence of the bare couplings for models 
with different values of A^ (defined keeping physical quantities such as the string 
tension fixed) becomes more accurate if one uses the mean field improved coupling 
(?mf proposed in Refs. ^ and obtained by dividing the lattice coupling by the 
average value of the plaquette, i.e. 



9l 



(A.2) 



Here we also consider the effective coupling (^cactus obtained within an all-order re- 



summation of cactus- type diagrams in perturbation theory and defined by 



2 

5'cactus 



9l 



(A.3) 



The function w{go) can be extracted by an appropriate algebraic equation that can 
be easily solved numerically: 



u e 



-u{N-l)/{2N) 



N-1 

N 



U 



(A.4) 



u{go) 



9l 



4(l-w(^7o)) 



where are the Laguerre polynomials. Data for A^ = 3, 4, 6 are shown in Fig. |^ 
(the data for A^ = 3 and some of those for A^ = 4 have been taken from Ref. |[49|| ) : 
we plot 1/7 = Ngl, l/7mf = A^5'mf and l/7cactus = Ngl^^^^ versus ad^/^ where a 
is the string tension. The gray points, corresponding to the mean field improved 
coupling, fall approximately on a single curve. A similar behaviour is observed for 
the cactus improved coupling (black points), while the data corresponding to the bare 
coupling (white points) show a wider spread as A^ is varied. We conclude that both 
the improved couplings provide an efficient tool to match theories at different values 
of A^. In this respect, the mean field improved coupling performs slightly better, 
but the advantage of the cactus definition is that it is determined from Eq. (|A.4|) , 
without requiring any simulation. 
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Figure 7: I/7, l/7mf and l/7cactus 
theories. 



versus oa^/^ fo^. SU(3), SU(4), and SU(6) lattice gauge 



B. Bulk phase transitions at large 

For sufficiently large values of and in particular in the large- limit, the Wilson 
lattice formulation of SU(A^) gauge theories presents a ffist order phase transition. 
This has been argued using various approaches, such as Monte Carlo simulations 
?5 , mean field calculations PUI , reduced models . In the following 



we present evidence for a ffist order transition at a finite bare coupling in the case 
of the SU(6) lattice gauge theory. On the other hand, in the SU(4) case no evidence 
of a bulk phase transition is found. 

B.l A first order transition for the SU(6) lattice gauge theory 

In the case = 6, we performed simulations starting from hot and cold configura- 
tions. Fig. ^ reports data for the energy density E (normalized to one for 7 = 0) 
obtained performing two simulations on a 8^ lattice, one starting from a hot config- 
uration and the other from a cold configuration. More precisely, in the ffist case we 
started from a hot configuration, and performed simulations starting from 7 = 0.334, 
increasing the value of 7 every 2000 heat bath updatings by 0.001. The data reported 
in the figure are the values of the energy obtained averaging over the second 1000 
updatings for each 7 value. The data concerning the simulation starting from a cold 
configuration were obtained similarly, i.e. starting from 7 = 0.344 and decreasing its 
value by 0.001 every 2000 upgradings. Fig. R| shows clearly the presence of hysteresis. 
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Figure 8: Energy density versus 7 obtained as explained in the text. The two vertical 
lines show the estimate of 7c obtained by the mixed-phase method. 

from which one may estimate 7c ~ 0.339. Note that the latent heat is relatively large 
indicating that the first order transition is rather strong. Another estimate of 7c can 
be obtained from the so-called mixed-phase method. One starts from a configuration 
that is half cold and half hot, and see which phase wins as 7 is varied. An estimate 
of 7c is obtained from the boundary values of 7 for which the change of final phase 
occurs. The two phases are easily recognized by their value of the energy: as shown 
in Fig. I, Ehot ^ 0.57 and E^oid ^ 0.47 at 7c. We obtained 7c ^ 0.3393, 7c ^ 0.3390 
and 7c ~ 0.3389 respectively from simulations on 8^, 10^ and 12*^ lattices. As a 
final estimate we consider 7c = 0.3389(4). For comparison, we mention the earlier 
estimate 7c = 0.333(14) 

B.2 Just a crossover for the SU(4) lattice gauge theory 

In the case of the SU(4) lattice gauge theory, no evidence for a bulk transition is 
observed. In particular, the specific heat does not appear to diverge with increasing 
lattice size. In Fig. ^ we plot data for the specific heat 

= (1 - ^(Tr f/,.(x))) (B.l) 

for various lattice sizes, i.e. 6^, 8^ and 12*^. The data were obtained from runs with 
typical statistics of 3-5x10*^ sweeps. They show a rather pronounced peak around 
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Figure 9: Specific heat versus 7 for 5C/(4) as obtained from simulations on 6^, 8^ and 
12^ lattices. Data for L = 6, 12 are slightly shifted horizontally to make the figure more 
readable. 



7 ~ 0.325 (corresponding to (3 = 10.4), but they appear to converge for increasing 
lattice size. We recall that, in the case of a first order phase transition, the finite-size 
scaling behavior of the peak value of the specific heat should diverge as 



C 



H,peak 



(B.2) 



thus Cn^peak ~ -^^ in our case. The data of Fig. |^ are clearly inconsistent with such 
a behavior. These results indicate that the Wilson formulation of the SU(4) lattice 
gauge theory does not undergo a first-order phase transition, but rather they suggest 
that, as in the SU(3) case, it exhibits a crossover between the strong and weak 
coupling regime, characterized by a rather pronounced peak of the specific heat. 



We mention that early works, based on mean field calculations and Monte 
Carlo simulations suggested a bulk first-order phase transition also for SU(4) 

and ~ 10.4. 



C. Critical slowing down for the topological modes 

Monte Carlo simulations of critical phenomena in statistical mechanics and of the 
continuum limit in quantum field theory are hampered by the problem of critical 
slowing down. For a general introduction to critical slowing down in Monte Carlo 
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simulations, see e.g. Ref. ||6^. The autocorrelation time r, which corresponds to 
the number of iterations needed to generate a new independent configuration, grows 
with increasing length scale ^; usually it blows up following a power law, i.e. r ~ 
In gauge theories one may consider = cr^^^^ as a length scale. Critical slowing 
down for traditional algorithms, such as standard Metropolis or heat bath, arises 
essentially from the fact that their updating is local: in a single step of the algo- 
rithm, information is transmitted from a given site/link to the nearest neighbors. 
Roughly, one expects that this information spreads following a random walk around 
the lattice. An essentially independent configuration is obtained when the informa- 
tion travels a distance of the order of the correlation length ^. This suggests that 
r ~ This guess is correct for Gaussian (free field) models; in general one expects 
that T ~ where 2; is a dynamical critical exponent. In the case of local algorithms, 
such as Metropolis and heat bath, one expects z ~ 2. Appropriate overrelaxation 
procedures may achieve a reduction of z, although the condition z > 1 holds for 
local algorithms. On the other hand, in the presence of relevant topological modes, 
the random-walk picture may fail. These modes may give rise to sizeable free-energy 
barriers separating different regions of the configuration space. The evolution in the 
configuration space may present a long-time relaxation due to transitions between 
different topological charge sectors. (Although lattice field theories cannot strictly 
possess topological properties, these are expected to be recovered in the continuum 



limit.) As in glass models, see e.g. Ref. and in liquids below a crossover transi- 



tion, see e.g. Ref. [Q, the presence of significant free-energy barriers may determine 
an effective separation of short-time relaxation within the free-energy basins from 
long-time relaxation related to the transitions between basins. The mechanism un- 
derlying this long-time relaxation is rather different from the simple random-walk 
spread of information, so the autocorrelation time of the topological modes may not 
show a simple power law behavior. This picture was also behind the so-called heating 
method |^ to measure the lattice renormalizations of topological charge operators. 



This method exploits the critical slowing down of the physical topological modes 
in off-equilibrium simulations, to disentangle them from the short-distance lattice 
renormalizat ions . 

Let us consider an observable O; its autocorrelation function Co{t) {t is the 
Monte Carlo time, i.e. the integer counting Monte Carlo iterations at equilibrium) 
is defined as 

Co{t) = {{0{t)-{0)){0{0)-{0))), (C.l) 

where the averages are taken at equilibrium. The integrated autocorrelation time tq 
associated with O is given by 
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Estimates of tq can be obtained by a blocking analysis of the data, without measuring 
the autocorrelation function Co{t). Indeed the following relation holds 

^o = iV^^, (C.3) 

where is the number of sweeps between two measurements of the observable O, 
Eq is the naive error calculated without taking into account the autocorrelations, and 
E is the correct error found after the blocking procedure (the estimate is meaningful 
only if A^in < To). Of course tq depends on the observable O; the largest value tq 
among the various observables provides the time scale to obtain a new independent 
configuration in simulations at equilibrium. In Ref. the autocorrelation time r^, 
of small Wilson loops was found to behave approximately as Ty^ ~ for the SU(3) 



gauge theories, see also Ref. As we shall see, a more severe form of critical 

slowing down is observed in measuring quantities related to the topological modes, 
such as the topological charge Q. 

We used cooling to determine Q from each lattice configuration, measuring Q 
every A^m = 100 Monte Carlo sweeps. Estimates of the (integrated) autocorrelation 
time Tq were obtained by a blocking analysis of the data, using Eq. (|C.3|) . In the 
SU(6) case we found tq ^ 268 for 7 = 0.342, tq ^ 466 for 7 = 0.344, tq ^ 1750 
for 7 = 0.348, and tq > 3000 for 7 = 0.350. Moreover we found that the run 
at the largest value of 7 considered, i.e. 7 = 0.354, did not correctly sample the 
topological charge presumably because the expected tq is large and the run was not 
sufficiently long (approximately 300k sweeps). In the SU(4) case we found tq ^ 101 
for 7 = 0.338, Tq ^ 210 for 7 = 0.341, and tq ^ 410,434 for 7 = 0.344 and 
L = 16, 12. The uncertainty on the above numbers should be at most 5% for SU(4) 
and 10-20% for SU(6). 

These estimates of tq are suggestive of an interesting phenomenon. As shown 
in Fig. |10|, they are consistent with an exponential critical slowing down. Indeed the 
data for the autocorrelation time can be well fitted by an exponential behavior of 
the type tq oc exp(c^o-) with c ~ 2.5 for SU(6) and c ~ 1.7 for SU{4:). One might 
also guess an increase of the constant c according to c ~ A^. We also mention that 
the data for tq are definitely inconsistent with a behavior of the type with z ~ 2. 
Indeed, an acceptable power law behavior fitting reasonably well the data for tq 
would require 2; ~ 7 for SU(4) and 2; ~ 9 for SU(6). We expect that a similar critical 
slowing down phenomenon occurs also for SU(3), and more generally in the presence 
of dynamical fermions. 

This dramatic effect was not observed for the other quantities considered in our 
study, such as the correlations determining the fc-string tensions and the glueball 
masses. A blocking analysis of the data did not show significant time correlations in 
measurements taken every 10-20 sweeps for all values of 7 considered. For instance, 
in the case of SU(4) and for 7 = 0.344 and L = 16 the autocorrelation time Tp 
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Figure 10: Plot of the autocorrelation time tq of the topological charge versus a 



-1/2 



of smeared and blocked Polyakov line correlators was estimated to be smaller than 
10, more precicely rp ^ 6, by a blocking analysis of the data. This fact suggests 
an approximate decoupling between the topological and nontopological modes. This 
seems to be also supported by the fact that string tension results for 7 = 0.354 in the 
SU(6) case (see Table extracted from a simulation which did not sample correctly 
Q, turn out to be in agreement with those for smaller 7, for which Q was sampled 
correctly. 

Such a critical slowing down phenomenon was already observed in simulations of 
two-dimensional CP'^"^ models |70[|. Given that these models possess some of the 
properties expected to hold in QCD |71], ^ (asymptotic freedom, a form of confine- 
ment due to the f/(l) gauge invariance, a nontrivial topological structure), they have 
been used as a laboratory to check and develop methods to investigate topological 
properties. Monte Carlo studies of lattice formulations of two-dimensional CP^~^ 
models, using local algorithms (mixtures of overrelaxation and heat bath upgradings, 
similar to the ones generally used for the four- dimensional SU(A^) gauge theories), 
have shown that the critical slowing down of the topological modes is consistent 
with an exponential growing with respect to the correlation length, worsening with 
increasing A^. On the other hand, nontopological quantities, such as the mass gap 
and Wilson loops, turned out not to be affected by this problem, suggesting a large 
decoupling between the topological and nontopological modes. 
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